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j^V We perform a study of the interfacial properties of a model suspension of hard sphere colloids with 

diameter o c and non-adsorbing ideal polymer coils with diameter o p . For the mixture in contact 
with a planar hard wall, we obtain from simulations the wall- fluid interfacial free energy, 7 w f , for size 
ratios q — a p /o c = 0.6 and f , using thermodynamic integration, and study the (excess) adsorption 
of colloids, r c , and of polymers, T p , at the hard wall. The interfacial tension of the free liquid-gas 
interface, is obtained following three different routes in simulations: i) from studying the system 
size dependence of the interfacial width according to the predictions of capillary wave theory, ii) from 
the probability distribution of the colloid density at coexistence in the grand canonical ensemble, 
and iii) for statepoints where the colloidal liquid wets the wall completely, from Young's equation 
55 ' relating 7i g to the difference of wall-liquid and wall-gas interfacial tensions, 7 w i — 7 wg . In addition, 

we calculate 7 w f , F c , and T p using density functional theory and a scaled particle theory based on 
free volume theory. Good agreement is found between the simulation results and those from density 
functional theory, while the results from scaled particle theory quantitatively deviate but reproduce 
some essential features. Simulation results for 7i g obtained from the three different routes are all 
, in good agreement. Density functional theory predicts 7i g with good accuracy for high polymer 

■ reservoir packing fractions, but yields deviations from the simulation results close to the critical 

O ' point. 
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I. INTRODUCTION 

Mixtures of sterically-stabilized colloids and non-adsorbing polymers are widely studied complex fluids Q, H ■ 
Provided the size and the number of polymers is sufficiently high, such mixtures can phase-separate into a colloidal 
gas phase that is poor in colloids and rich in polymers, and a colloidal liquid phase that is rich in colloids and 
poor in polymers. The mechanism behind this demixing transition is of entropic origin and is due to the so-called 
depletion effect: An effective attraction between the colloids is generated due to an unbalanced osmotic pressure 
arising from exclusion of polymer coils in the depletion zones between the colloids 0, ||| . Since the associated relevant 
i-^J ■ time and length scales are much larger than in atomic and molecular systems, direct experimental observations using 
advanced microscopy techniques enable the study of many interesting physical phenomena, e.g. thermal capillary 
O ■ waves at fluid interfaces and droplet coalescence were observed recently in real space and real time using confocal 
^ ' microscopy Other recent examples are the direct measurement of the contact angle of the colloidal gas-liquid 
J> , interface and different substrates , complete wetting of a substrate [p, || , a wetting transition from complete to 
partial wetting , and capillary condensation in confining geometry |lOj . Hence it is fair to say that these mixtures 
serve as excellent model systems to investigate fundamental issues in statistical physics. A particularly simple model 
for colloid-polymer mixtures was proposed by Vrij p|, and is often referred to as the Asakura-Oosawa-Vrij (AOV) 
model. A good historical introduction to this model together with many recent results can be found in the paper by 
Brader et al. la] . The bulk phase behavior as well as inhomogeneous properties of the AOV model were studied with 
both theory GH 111 III 111 and computer simulation [H TTrl 111 111 fill El 151 Ej |2^| . Recently, attention was 
paid to more realistic model interactions [H ll8L |24| for colloid-polymer mixtures. However, most of the essential 
physics of real colloid-polymer mixtures is indeed captured by the AOV model. 
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In this work we study the wall-fluid and liquid-gas interfacial tensions of the AOV colloid-polymer mixture using 
Monte Carlo simulations and check the predictions of density functional theory (DFT) |14L l2q based on an extension 
of the Rosenfeld functional pfij . and of a scaled particle theory (SPT) 0| based on the free volume theory [Tlj . 
The wall-fluid interfacial tension of a hard-sphere fluid in contact with a planar hard wall was recently calculated by 
Heni and Lowen using a thermodynamic integration procedure along a path that corresponds to the growth of a wall 
in a bulk system [2?1 |. Here, we propose a thermodynamic integration approach similar in spirit, to determine the 
wall-fluid interfacial free energy of the AOV colloid-polymer mixture in contact with a planar hard wall. 

Different routes exist to determine the liquid-gas interfacial tension. The pressure tensor metho d l28t |29| is par- 
ticularly suitable for Molecular Dynamics simulation studies. The probability distribution method [3(j was used in 
Monte-Carlo simulations ElO^l and was recently applied to calculate the liquid-gas interfacial tension of colloid- 
polymer mixtures pH l22l l23t l33 | . The authors showed that the interfacial width depends on system size and they 
verified the predictions from capillary wave theory on system size dependence |22| . Sides et al. and Lacasse et al. found 
good agreement between the interfacial tension obtained from the pressure tensor method and from the predictions 
of capillary wave theory provided that the density profile is fit to an error function for Lennard-Jones particles and 
polymer blends l35l ] . In this work, we use both the capillary wave theory, as proposed in Refs. |34l Is^l] . as well 
as the probability distribution method to calculate the liquid-gas interfacial tension. In addition, we employ Young's 
equation in the complete wetting regime to obtain an estimate for the tension of the liquid-gas interface from the 
difference in wall tensions, and compare our results with DFT calculations. 

The paper is organized as follows. In section [H] we briefly review the AOV model. In section UlI Al an overview of 
the thermodynamics of inhomogeneous systems is given. In section IIII Bl we lay out the simulation methods used in 
the determination of the wall- fluid tension. Section fill CI details the simulation method to obtain the adsorption at a 
hard wall. Section [ill Dl presents the corresponding derivation using scaled particle theory. In section Till El we review 
the capillary wave theory for the width of a liquid-gas interface. Sections IIII Fl and IIII CI are devoted to obtaining the 
liquid-gas interfacial tension from the probability distribution of the colloid density at coexistence and from Young's 
equation, respectively. Section Till HI gives a brief overview of the DFT. In section ITVl we discuss our results for the 
wall- fluid tension (|I V Al) . adsorption (|IV Bfl . and liquid- gas interfacial tension IjIV Cfl . Concluding remarks are given 
in section Ivl 



II. MODEL 



We consider a mixture of sterically-stabilized colloidal spheres (species c) and non-adsorbing polymer coils (species 
p) immersed in a solvent. The interaction between two sterically stabilized colloids resembles closely that of hard 
spheres, while a dilute solution of polymer coils in a theta-solvent can be treated as an ideal gas as the polymer 
coils are interpenetrable and noninteracting. The polymer coils are assumed to be excluded from the colloids to a 
center-of-mass distance of (er c + d p )/2, where a c is the diameter of the colloids, and a p = 2R g is the diameter of the 
polymer coils, with R g being the radius of gyration of the polymers. A simple idealized model for such a mixture is 
the so-called Asakura-Oosawa-Vrij (AOV) model 0, Q defined through pair potentials, that between colloids being 

m \ _ / 00 if Rij < °c m 

vecyttij) - <y otherwisej W 

where Rij = \Ri — Rj\ is the center-of-mass distance between two colloidal particles, with Ri (Rj) the center-of-mass 
coordinate of colloid i (j). The polymers are described as noninteracting particles with 

Vppinj) = 0, (2) 

where = \r*i — fj\ is the distance between two polymers, with f*j (fj) the center-of-mass coordinate of polymer i 
(j). The colloid-polymer interaction is 

v04-r iD -{- "MA^I <<*+*>/» (3) 

where \Ri — fj\ is the distance between colloid i and polymer j. Our simulations are performed in a box with 
dimensions L x L x H with three-dimensional periodic boundary conditions in the case of bulk simulations. To create 
a wall-fluid or a liquid-gas interface, we perform simulations in a box with periodic boundary condition solely in the 
x and y directions and two impenetrable hard walls in the z direction. The wall-particle potential acting on particles 
of species k = c,p reads 

, , ( if cr fc /2 < z Kl <H k - cr fe /2 , . 

v w k{z k ,i) - s ^ otherwise , KV 
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where Zf.,% is the z-coordinate of particle i of species k, and Hu is the separation distance between the two walls for 
species k. For the simulations of hard wall properties, we use H c = H p = H, corresponding to two planar hard 
walls. For the simulations of the liquid-gas interface we use a box with periodic boundary conditions in the x and y 
directions and delimited in the z direction by one impenetrable and one semi-permeable wall; this asymmetric slit is 
defined by the wall-particle potential |@} with H c = H — 2a c and H p = H . The impenetrable wall at z = favors 
the colloidal liquid phase, while the semi-permeable wall at z — H c is impenetrable for the colloidal particles, but 
penetrable for the polymers (which are free to overlap with this wall). Hence, there is no effective polymer-mediated 
wall-colloid attraction and the semi-permeable wall favors the colloidal gas phase 0, 0, 0, |2(j ■ 

The size ratio q = o p ja c is a geometric parameter that controls the range of the effective depletion interaction 
between the colloids. Packing fraction are denoted by % = (tto-IN^)/ (6L 2 H), and the number density is denoted by 
Pk = Nk/(L 2 H) for species k = c,p. We also employ, as an alternative thermodynamic variable to rj p , the polymer 
reservoir packing fraction 

Vp = Q^pPp = g°p^ (5) 
where z p is the fugacity of the polymers, that constitute an ideal gas with density p r in the reservoir. 



III. METHODS 



A. Overview of interfacial thermodynamics 



Generally, the interfacial tension in an inhomogeneous system is the grand potential per unit area needed to create 
an interface in an initially uniform bulk system at fixed chemical potential of colloids, /i c , and polymers, ji p , and fixed 
volume V and temperature T. The grand potential for a bulk mixture of colloids and polymers reads 

ft bul V, fx p , V, T) = -p(n c , p p , T)V, (6) 

where p is the bulk pressure. The system in contact with an interface possesses the grand potential 

0(/x c , fip, V, T, A) = ~p(p c , p p , T)V + 7(/i c , /J, p , T)A, (7) 

where A is the area of the interface and 7(/i c , T) is the interfacial tension, which can hence be expressed as 

nfrg, p p , v, t, A) - n bulk ( Pc , tip, v, t) 



A 



(8) 



Besides the liquid-gas interface, where 7 = 7i g , Eqs. ©, 0, and JBJl apply also for a fluid adsorbed between two 
parallel plates (walls), where 7 = 7^, provided that the wall separation is sufficiently large j3||[37l, and that the area 
A is equal to the total area of the two plates, A — 21? . At fixed chemical potentials the number of particles in the 
inhomogeneous system, N c and N p , of colloids and polymers, respectively, will be in general different from those in 
the bulk, N bulk and N p ulk . The excess number of colloids and polymers per unit area, i.e. the adsorptions T c and 
T p , respectively, are denned as 

N c — A bulk 

T c (p c ,p p ,T) = — , (9) 

N p - A bulk 

r p (p c ,p p ,T) = A p . (10) 

The grand potentials 10 and Q) in differential form read 

dn h » lk (p c ,p p: V 7 T) = -A c bulk d Mc - N™ k dfip-pdV-S hu)k dT,, (11) 
dn{p c ,p p ,V,T,A) = -N c dp c - N p dp p -pdV - SdT + -fdA. (12) 

Using Eqs. I|ll|) and l|12l) and Eq. (JSJl in differential form, it is straightforward to show (3^ that the adsorptions are 
related to the interfacial tension through 

r* = -(P-) and T P = -(£-) . (13) 
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B. Hard wall-fluid interfacial tension via thermodynamic integration 



To determine, from simulations, the wall-fluid tension 7 w f of the AOV model we should apply equation ijSJl. as is 
manifest in the grand canonical ensemble, i.e. for constant colloid and polymer fugacities. However, in our simulation 
it is more convenient to use the semi-grand canonical ensemble fixing the number of colloids and the fugacity of the 
polymers. The reason is twofold. First the interfacial tension as a function of the fugacity of polymers can be directly 
compared to the DFT results of Ref. 01 • Second, fixing the number of colloids instead of their fugacity allow us to 
efficiently study state points with high packing fractions of colloids; generally grand ensemble simulations are difficult 
to perform at high densities due to small particle insertion probabilities. To compute the tension we have to recast 
Eq. JSJ in a way that is consistent with the semi-grand canonical ensemble. The grand potentials for the bulk and 
the inhomogeneous system are related to the corresponding Hclmholtz free energies via a Legendre transformation, 

O bul V, fi p , V, T) = F bulk (AT bulk , iV bulk , V, T) - Mc 7V bulk - fi p N* ulk , (14) 
n(/i c ,/x P ,V,T,A) = F(N C ,N P ,V,T,A) - n c N c - v P N p . (15) 

We substitute Eqs. (fTl|) and (fTSf) in Eq. JSJl to obtain 

F(N c ,N p )-F b ^ k (N^ k ,N^ k ) 

i = 7 M c r c - Ai P r p , (16) 

where we omitted the dependence on the variables V, T, \i c , and \x p in the notation. Note that the tension is not only 
the difference of the Helmholtz free energies, but additional terms, /x c r c and n p T p , arise in Eq. l|ltj|) . One can further 
simplify by Taylor expanding F(N C , N p , V, T, A) around N^ ulk : 

BF 

F(N C , N p , V, T, A) = F(N* , N p , V, T, A) + — (N c - 7V bulk ) + 0{{N C - iV bulk ) 2 ). (17) 

Keeping only the first order term, one can approximate the interfacial tension as 

F(iV bulk , N p , V, T, A) - F hulk (N^ ulk , AT bulk , V, T) 
7 - j M P r p . (18) 

The same approximation was employed in Ref. [27j (using N p = and iVp Ulk = 0) to calculate the interfacial free 
energy of hard spheres in contact with a planar hard wall. To compute the wall tension, we need to perform two free 
energy calculations, one for the bulk and one for the inhomogeneous system. As the free energy cannot be measured 
directly in a Monte Carlo simulation, we use the thermodynamic integration technique |39|| to relate the free energy 
of the system of interest to that of a reference system 

F(N C , N p , V, T,A,X= A max ) = F ld (N c , N p , V, T, A = 0) + f ™" dx(^-) . (19) 



A=0 



The reference system is chosen to be an ideal gas, so Fid(N c , N p , V, T, X = 0) is the Hclmholtz free energy of N c ideal 
colloids and N p ideal polymers in a volume V at temperature T . We then introduce the suitable auxiliary Hamiltonian 

N c N c N p N c N p \ 

H X = X\^2 V cc (Rij) Vc p^ - F iD +^ V ^c{z c ,i) + e£ V„ P (2p,i) , (20) 

^i<j i=l j=l i=l i=l J 

where we approximate the hard-core potentials of the AOV model with penetrable potentials: The colloid-colloid 
interaction reads 

V cc {Rij) = e(<r c - IU-), (21) 

where 0(x) is the Heaviside step function. Likewise we define the interaction potential between the colloids and the 
polymers as 

V cp {\Ik - r 3 \) = 9 (^l±^l - IR. _ ^.|) . (22) 
The interaction between the walls and the particles of species k = c, p is 

V wk (z k>i ) = © (y - «fc,<) + Q ( y - (-fffc - z u)) , (23) 
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where Zk,% is the z-coordinate of particle i of species k. For A = 0, this system reduces to an ideal gas, while for 
A — * oo, the system describes the AOV model given by Eqs. Q~0 in bulk (e = 0) or confined by two walls (e = 1). 
The interaction potential is switched on adiabatically using the coupling parameter A. In principle, our system of 
interest is described by the Hamiltonian (|20[) using A max — > oo, but also for sufficiently high values of A max the system 
reduces to our system of interest with hard-core potentials. Clearly, A max should be sufficiently large to ensure that 
the system is indeed behaving as the hard-core system of interest. On the other hand, A max should not be too large, 
as this would make the numerical integration less accurate. The integrand function of Eq. i|19|l is then 

The function (^) x is computed counting the number of overlaps between colloids, colloids and polymers and (for the 
inhomogeneous system only) particles and walls. The free energy can then be obtained using Eq. The integrals 

are evaluated using a 21-point Gauss-Kronrod formula, where 5000-15000 MC cycles per particle are used for the 
sampling of each integration point. The wall- fluid interfacial tension is then computed using equation (TSJ. In detail, 
we first perform a simulation of the AOV model in bulk and a separate simulation of the AOV model confined by two 
walls, given by the interactions 0J-J2J, both in the semi-grand canonical ensemble, i.e., we fix the number of colloids 
N~ , the chemical potential of the polymers fi p , and the volume V = L x L x H . We measure the average number 
of polymer in the bulk, (A^ ulk ), and in the confined system, (N p ). We then perform two separate thermodynamic 
integrations (in the canonical ensemble) to obtain the free energy of the bulk system with A^ ulk colloids and (iVp ulk ) 
polymers in a volume V, and the confined system of volume V with ]\[^ ulk colloids and (N p ) polymers. In the canonical 
ensemble simulations, we determined the chemical potential of the polymer as a consistency check. Typical number 
of the colloids and the polymers are N]? = 54 — 900 and ulk = — 20000, while the volume of the simulation box 
is about V = (1200 — 3000)erj? and H > 16cr c . The errors are estimated calculating the standard deviation from 4 or 
5 independent simulations. 



C. Adsorption at a hard wall from simulation 



To study the (excess) adsorption of colloid-polymer mixtures at a planar hard wall we simulated both the bulk 
mixture and the mixture in contact with the hard wall in two independent Monte Carlo simulations in the grand 
canonical ensemble ensemble and hence we considered only statepoints of low colloid packing fraction r\ c . After 
discarding 50000 MC steps per particle for equilibration, we take the average of the number of particles for another 
50000 MC steps per particle. The differences in particle numbers (per unit area) in the confined system and in the 
bulk system then give the adsorption of both species via Eqs. © and l|10|) . 



D. Adsorption at a hard wall from scaled-particle theory 



For a system of hard spheres the scaled particle theory |40j, |41| describes quite accurately the pressure p, the hard 
wall-fluid interfacial tension 7h s , and the (excess) adsorption Ths, given through the expressions 

0p _ l + ri c + ri 2 c 

lhsCT c = - — r . (27) 

7r(l + 2t? c ) 7T 

In particular, Eq. I|27|l was shown to compare well with simulation |2^| and DFT (4^ results. 

Recently an SPT expression for the wall-fluid tension of AOV model colloid-polymer mixtures was derived by 
Wessels et al. llj using the bulk free energy for a ternary mixture obtained from free volume theory as an input, 
and taking the limit of vanishing concentration and infinite size of the third component. Their expression reads 

/?7wf<7 c 2 = /? 7 h s ^ 2 +^/(r?c), (28) 

where f(rj c ) = 3a(?7 c )/(g 2 7r)[l + (1 + 3q + q 2 )r + (3q + 4g 2 )r 2 + 3q 2 T 3 ], r = T) c /(1 — r] c ) and /37h s c 2 is given by equation 
l|26p. The polymer free volume is given by the scaled particle theory as a(rj e ) = (1 — r/ c ) cxp( — (3q + 3q 2 + q 3 )r — 
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(9g 2 /2 + 3<7 3 )t 2 — 3g 3 r 3 ). Results for 7 w f from equation l|28[l were found in Ref. 01 to compare reasonably well 
with those from full numerical density functional calculations. Below we will compare these approaches against our 
simulation data. 

Here we derive an SPT expression for the adsorption of the AOV model at a hard wall starting from equation l|28l) 
and building derivatives according to (| 1 31) . The colloid chemical potential obtained from the free volume theory [llj 
is 



(3 fl c = PflhsiVc) ~Vp-^, 



(29) 



where PfihsiVc) = ?7c(14 — 13r/ c + 57y 2 )/(2(l — r/ c ) 3 ) — log(l — r/ c ) + log(6r/ c /7r) is the SPT expression of the chemical 
potential of a system of pure hard spheres at packing fraction r/ c and a = da/drj c . We compute the colloidal 
adsorption using equation 113fl 



9 / 97wf(^c,Mp) cr c _ df3"f w {(rj c ,^ p )(j 2 c dr\ c 



where 



d{3[i c drj, 



d(3 He 



is computed using equation l|29|) . The final expression is 



dr] c 



1 I jf—L — \ 
1~% 



(30) 



(31) 



(32) 



where / = df/drj c , 7 hs = djh s /di] c , fi hs = d/j,hs/dr) c and a = d 2 a/drj 2 . We note that the hard-sphere limit is 
obtained correctly for 77^ = 0. We also calculate the polymer adsorption 



r 2 d(3-/ w {(r] c ,ii p )a 2 



= -v r J 



drjc 



r df(n 
'ip 



drj c 



di] c J dfi[ip 



Rewriting Eq. l(2*9"|) as 



V P 



PfJ-hsjVc) ~ PfJ-c 3 



(33) 



(34) 



we arrive at 



The final expression reads 



drj c 



drip 



r l 

rj p a 



tf 3 A»L - % a " 



(35) 



r p(Vc)Vc = -(^7wf - Plhs)cr 2 c + rf p —T c (r] c )a 2 c 



(36) 



E. Liquid-gas interfacial tension from capillary wave broadening 



Capillary wave theory describes the broadening of an intrinsic interface of width wo due to thermal fluctuations. 
This broadening depends primarily on the interfacial tension and the area of the interface. To calculate the capillary 
wave contribution to the interfacial width one has to sum over the contributions from each thermally excited capillary 
wave to the amplitude of the oscillations in the instantaneous interface position. Here we briefly sketch the derivation 
of Sides et al. and Lacasse et al. and refer the reader to Refs. [34ll35| for further details; a very recent study devoted 
to capillary waves in colloid-polymer mixtures is that by Vink et al. |43j |. Fluctuations due to capillary waves in 
C(x,y), the mean location of the interface in the z-direction, have an energy cost due to the increase in surface area 
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of the interface. The free energy cost of the interfacial fluctuations is the product of the excess area of the undulated 
interface over that of the flat one, and a liquid-gas interfacial tension 7i g , which is assumed to be independent of 
curvature. The interfacial Hamiltonian, assuming that £ and its derivatives are small, reads 

n = ^ [ dxdy\V((x,y)\ 2 . (37) 



2 

Introducing a Fourier expansion of C(%, y), one arrives at 

n = "2T / d ^ 2 l^l 2 ' ( 38 ) 

where q = (q x ,q y ) denotes a two-dimensional wavevector and £(<f) is the Fourier transform of Using the 

equipartition theorem, the mean-square amplitude for each interfacial excitation mode reads 

The mean-squared real space fluctuations can be calculated by summing over all allowed modes: 



1 

(2*) 

k B T . / L 



(IC(^)I 2 ) = 7^2 / dqdCm 2 } 



2 V n UJ' (40) 

where the low q cutoff, q m i n , is determined by the system size, i.e., g m ; n = 2~k / L in our simulations and gives rise to 
system size dependence. The high q cutoff, q max — 27r/^, is determined by the bulk correlation length £t>, which is of 
the order of the colloid diameter, and avoids the divergence in the integral. 

The total width of the interface, as measured in experiments and simulations, includes contributions from the intrin- 
sic width and the broadening due to capillary wave fluctuations. If one assumes that the capillary-wave fluctuations 
are decoupled from the intrinsic profile, the total interfacial profile *&(z) can be expressed as a convolution of the 
intrinsic interfacial profile ip{ z ) an d the fluctuations due to capillary waves 

/•OO 

tf(z)= / - z )P{z )dz , (41) 



where P(zq) is the probability of finding the local interface at zq. The interfacial order parameter profile ^f(z) is 
defined such that it varies between -1 and 1 



*(*) = "I 9 

Pc ~ Pc 



f \ P\ + P 9 c 
Pc{z) S 



(42) 



where p c (z) is the cross-section averaged density profile of the colloids and p\ and p 9 c are the colloid densities of the 
colloidal "liquid" and "gas" phase at coexistence. In Ref. 0,113 > the authors define the variance of the derivative of 
the total interfacial profile d^(z)/dz = ^' as a measure of the width of the interface. The variance of a distribution 
/ reads 

_ Sj^fmdz -(dydq*)f(q)\ q =o 

v[n - JZcMdz - /(o) ' { 6) 

where f{q) is the Fourier transform of f(z). Using this choice of measure for the interfacial width, one can show 
explicitly using the convolution theorem that the total interfacial width can be written as the sum of an intrinsic part 
and a contribution due to capillary wave fluctuations 

v[V'] = v[ip'] +v[P] 

k B T , ( L 



= *'i + 2^ 1 H& 1 ' (44) 



where we identify v[P] as the mean-squared fluctuations due to capillary waves, i.e., (\((x,y)\ 2 ). 



8 



To measure the tension using the results of capillary wave theory one needs to create a liquid-gas interface in the 
simulation box. To stabilize the liquid-gas interface we perform the simulations in a box with dimensions L x L x H, 
with H = 48er c , delimited in the z direction by one impenetrable wall and one semi-permeable wall. We vary the 
lateral dimensions in the range of 5er c < L < 25er c . The canonical simulations are started in the middle of the two- 
phase region. After equilibration the system is phase separated, the gas phase in contact with the semi-permeable 
wall and the liquid phase in contact with the hard-wall. The liquid-gas interface is in the middle of the simulation 
box (see also Fig. [HJl. We determine the total interfacial width v[*f?'] from the interfacial order parameter profile ^f(z) 
measured in simulations according to Eq. I|42|l . A priori, it is not clear whether the density profile should be fit to an 
error function or a hyperbolic tangent. In Ref. j3^|, it is shown explicitly that the interfacial width extracted through 
fits of the density profile to a hyperbolic tangent lead to systematic errors and that better results for the interfacial 
tensions are obtained using the error function. In addition, later work on water/carbon tetrachloride |44| and molten 
salt (KI) 0| interfaces also found good agreement between the interfacial tension calculated using the capillary wave 
formalism and the pressure tensor components. Fitting ^(z) by an error function eri[(z — zq)/ (wy/2)], using zq and 
w as fitting parameters, we find that the variance of the derivative of this fitting function is related to the interfacial 
width w, i.e., v [*'] = w 2 . Using Eq. JH we are able to determine j\ g from the fits of the size dependence of the 
interfacial width. 

In contrast to the wall- fluid tension and the adsorption simulations described in section IIII Bl and IIII CI we used 
for the liquid-gas interfacial tension simulations described in this section and in section IIII Fl an efficient simulation 
scheme for the AOV model that was recently developed ^ is based on the exact effective one-component 

Hamiltonian of the colloids, i.e., it incorporates all the effective polymer-mediated many-body interactions. The 
effective one-component Hamiltonian can be derived by integrating out the polymer degrees of freedom in the binary 
colloid-polymer mixture. To this end, we consider N c colloids and N p polymer coils in a macroscopic volume V at 
temperature T. The total Hamiltonian consists of interaction terms H — H cc + H cp + H pp: where H cc = Yl^Jj v cc (Rij ), 

H C p = Y^h " 12^ P u cp(|Rj — r j | ) , and H pp — Yl!i<j v pp{ r ij)- It is convenient to consider the system in the (N c , z p , V, T) 
ensemble, in which the polymer fugacity is fixed. The thermodynamic potential F of this system can be written 
as exp[— j3F\ = Tr c exp[— f3H e ff], where H e ft = H cc + 17 is the effective Hamiltonian and where Tr c is short for 
l/JV c !Ag times the volume integral over the coordinates of the colloids, and where A c is the thermal wavelength. It 
is straightforward to show that one obtains for the present model the exact result /3f2 = —z p Vf, with the so-called free 
volume Vf, i.e., the accessible volume for the center-of-mass of the polymer coils. This free volume can be calculated 
numericall y o n a smart grid for each static colloid configuration. For more details, we like to refer the reader to 
Ref. |17l |49 | . The advantage of this scheme is that the polymer degrees of freedom are integrated out and enter 
the effective one-component colloid Hamiltonian only by the polymer reservoir packing fraction rf . Hence, we avoid 
equilibration and statistical accuracy problems due to fluctuating polymer numbers. Moreover, our simulations are 
not limited by the total number of polymers and can be performed at high polymer reservoir packing fractions far 
away from the critical point. 

F. Liquid-gas interfacial tension from the probability distribution 

In addition, we determine 7i g using the probability distribution method. We perform Monte-Carlo simulations in 
the grand canonical ensemble using again our novel effective one-component simulation scheme [T^ |. explained in Sec. 
IIII El . To obtain the probability P(N c )\ Zctri r of observing N c colloids in a volume V at fixed colloid fugacity z c and 
fixed polymer reservoir packing fraction jjl, we use a sampling technique called successive umbrella sampling |47| . 
Employing this technique, we sample the probability distribution P(AT c )| Xo „r in small windows one after the other, 
in which the number of colloids N c is allowed to fluctuate between and 1 in the first window, 1 and 2 in the second 
window, etc. We first perform an exploratory short run without a bias, yielding P{N C — 1) and P(N C ) for window 
N c . We then perform a biased simulation in which we sample from a non-physical distribution ir(N c ) = g(N c )P(N c ) 
instead of the grand canonical distribution P(N C ). We have chosen the weight function g{N c ) = 1/P(N C ), where 
we use P(N C ) obtained from the unbiased exploratory simulation run. This choice for the weight function yields a 
constant constant ir(N c ) and the system will visit equally the state with N c — 1 colloids as the state with N c colloids. 
Of course, we have to correct our weighted sampling for the bias by dividing out the weight function g(N c ). The more 
accurate grand canonical distribution P(N C ) is obtained from 



(45) 
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In addition, we use the histogram reweighting technique to obtain the probability distribution for any z' c once 
P(N c )\ Zc<ri r is known for a given z c |33| : 

P(Nc)k K = lnP(N c )\ Zc , % + In N c . (46) 

At phase coexistence, the distribution function P(N C ) becomes bimodal, with two separate peaks of equal area for 
the "colloidal " liquid and gas phase. To determine phase coexistence, we normalize P(N c )\ ZctV r to unity 



P(N c )\ Zc , % dN c = 1 (47) 
and we determine the average number of colloids 

(N c ) = / N c P{N c )\ Zc , vl dN c . (48) 



Using the histogram reweighting technique l|46(l we determine for which z' c the equal area rule 

-(N e ) 



P(N c )\ z ,, % dN c = P(N c )\ z , i% dN c , (49) 

J(N C ) 

representing the condition for phase coexistence, is satisfied. 

The liquid-gas interfacial tension 7i g for a finite system of volume V — L 2 H can be obtained from P(N C )\ Z / iV r at 
coexistence: 



^s,l = 7772 



1 

2L 2 



.max/ \ i/7-i/7vr \\ 

m I y ) ~~ ln \ P ( N c,min)) 



(50) 



where P{N-! max ) and P (-A^ max ) are the maxima of the gas and liquid peaks, respectively, and P(^V c , min ) is the 
minimum between the two peaks. We can determine the interfacial tension for the infinite system, i.e., 7i g , by 
performing simulations for a range of systems sizes and by extrapolating the results to the infinite system, as shown 
by Binder, using the relation [3(1 13lll33| 



xlnL InA 



where A and x are generally unknown. 



G. Liquid-gas interfacial tension from Young's equation 

Young's equation relates 7i g to the difference in wall-fluid tension for the gas, 7 wg , and the liquid phase, 7 w i, via 

(7w g - 7wi) = 7ig cc, s , (52) 

where 9 is the (macroscopic) contact angle at which the gas-liquid interface hits the wall. In the region of complete 
wetting, the contact angle is zero and hence 7i g can be obtained from the difference of the wall-gas and wall-liquid 
tensions, 7i g = (7 wg - j w \). 



H. Density functional theory for interfacial properties 

We use the approximation for the Helmholtz excess free energy for the AOV model as given in [25| . For given external 
potential, the density functional is numerically minimized using a standard iteration procedure. The interfacial tension 
is then obtained from Eq. ©, and the adsorption of both species from Eqs. 10, 1)11)0. Technical details about the 
DFT implementation that also apply to the present study are given in |l3j . 
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IV. RESULTS 
A. Wall-fluid interfacial tension 

In this section we present the results on the interfacial tension of model colloid-polymer mixtures from simulation, 
SPT and DFT. We checked the simulation technique performing simulations for a system of pure hard spheres 
(Vp = 0) an d started by finding the optimal value for A max , the maximum height of the potential l|20l) . To this end, 
we computed the average number of overlaps among particles and walls for different values of the potential 

height A, with H defined by the equation 120|) . The value of A max depends on the particle packing fraction. In Fig. 
^we plot the average number of overlaps between the colloidal particles and the walls for packing fraction rj c = 0.4. 
For all values of /3A > 10, with /3 = 1/kT the inverse temperature, the average number of overlaps is zero within the 
statistical fluctuations. Smaller packing fractions require smaller values of /3A max . We checked the reliability of the 
approximation by computing the reduced wall-fluid interfacial tension Pj^ga^. (Fig. |2J) and found that our simulation 
results are consistent with the results of Ref. |27|. At packing fraction r/ c > 0.45 the precision of the simulation is low 
but comparable with the simulation method of Heni and Lowen |2?| using the wall insertion method, and of Henderson 
and van Swol using the pressure tensor method psl ] . The agreement between simulation and density functional theory 
(DFT) 0112 is remarkably good. The scaled particle theory (SPT) Q overestimates the wall-fluid tension at high 
density. This is due to the the inaccuracy of the SPT equation of state. Better results can be obtained combining the 
SPT equation for the interfacial tension and the Carnahan-Starling equation of state as explained in detail in Ref. 

We now determine the wall-fluid interfacial tension for AOV colloid-polymer mixtures of size ratio q = 0.6 and 
q = 1 for different values of the polymer reservoir packing fraction 77^ and of the colloid packing fraction r\ c . The 
addition of nonadsorbing polymers to a colloidal suspension of hard-sphere can induce a phase separation. In Fig. 
121 we show the bulk phase diagram for size ratio q = 1 from previous simulations [l7| in the (rjp,rj c ) representation. 
For comparison, we also plot the phase diagram obtained from free volume theory, which is equivalent to our DFT 
phase diagram [l_l] . At rf p — we find the freezing transition of the pure hard-sphere system with packing fractions 

r\{ ~ 0.494 and rf c ~ 0.545 for the coexisting fluid and solid phase, respectively. The critical point is estimated to be 
at rfp Crit = 0.86, while DFT, equivalent to the free-volume theory predicts i]p Crit — 0.638. For rf v < eta r p crit , there is 
a stable fluid phase for r\ c < 0.494, a fluid-solid coexistence region for 0.494 < r\ c < 0.545, and a stable solid phase 
(fee crystal) for r\ c > 0.545. For rf v > rf crit , a fluid-fluid coexistence region appears where the system demixes in a 
colloidal liquid phase, rich in colloids and poor in polymers, and a colloidal gas phase, that is poor in colloids and rich 
in polymers. The triple point, where the gas, the liquid, and the solid are in coexistence, is located at ??p triplc = 6. 
For 7/p > i]p triple ^ ne fluid-fluid coexistence region disappears, and a wide crystal- fluid coexistence region appears. 
The overall phase diagram is analogous to that of a simple fluids upon identifying rf v with the inverse temperature. 
Despite differences near the critical point, DFT and simulations results are in good agreement for state points at 
> 1.5. In Fig. 0Ji and Fig. ^ we show the wall- fluid tension for state points below the gas-liquid critical point for 
size ratio q = 0.6 and q — 1, respectively. For comparison, we also plot the results for pure hard spheres (77^ = 0). The 
addition of non-adsorbing polymers to a suspension of hard-sphere colloids (i.e. increasing 77^) increases the wall-fluid 
interfacial tension. For r\ c = 0, the wall-tension is the work done to introduce an impenetrable wall in an ideal gas 
of polymers divided by the total area: pj(r] c = 0) = pP ld a p /2, where j3P td = p r p is the bulk pressure of the ideal 
gas of polymer and a p /2 is the thickness of the depletion layer of the polymer at the wall. For small r] c , the slope 
of the tension is smaller than in the hard sphere case and for 77^ >= 0.4 it is negative. This is due to the attractive 
interaction that arises between the colloidal particles and the walls. For large rj c the interfacial tension approaches 
that of pure hard spheres as at high colloid density the number of polymers in the mixture rapidly approaches zero. 
Simulations and DFT are in good agreement for all state points that we considered. The SPT predicts correctly the 
value at rj c — 0, but it systematically overestimates the wall-fluid tensions for all values of r\ c > 0. One can show that 
the low rjc expansion violates an exact wall sum rule |48| . The deviation increases with increasing 77^. In Fig. 0b we 
show the results for size ratio q = 0.6 for state points that are at higher 77^ than the DFT gas-liquid critical point. 
We did not calculate the binodal with computer simulation, but the system was still in the one phase region of the 
phase diagram for 77^ = 0.5 and 0.6. For comparison, we also plot the results from DFT and SPT. Note that DFT 
results are only shown in the stable gas and liquid regimes and are hence disconnected from each other, showing the 
biphasic region at intermediate ?7 C . In Fig.^Ji we show the results for the size ratio q = 1 for state points that are at 
higher 77^ than the DFT gas-liquid critical point. For small 77,, the SPT fails to reproduce the slope of the curves, due 
to the absence of colloid correlations (layering) near the hard-wall in SPT theory. 
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B. Adsorption at a hard wall 



In this section we present results on the adsorption of colloids and polymer at a hard-wall. We compare the 
simulation results with those from DFT and SPT Eqs. © and (|10fl . In Figs.[S^ and[Sp, we show the results on the 
colloidal adsorption while in Figs.JSJ: and[5Ji, we show the results on polymer adsorption for size ratio q=0.6 and q=l, 
respectively. We notice that increasing the number of polymer in the system (i.e. increasing ryp the adsorption of 
colloids increases; the colloids are attracted at the hard wall by the depletion interaction. As shown by the polymer 
adsorption the increase in number of colloidal particle at the walls is followed by a decrease of the number of adsorbed 
polymer while increasing the total number of polymers in the system. The agreement between simulations and DFT is 
good. This is not surprising since the DFT is known to provide an accurate description of the colloid-polymer mixture 
at a planar hard- wall j^j. The SPT equations reproduce the r] c =0 limit correctly. For r\ c ^ the essential features 
are reproduced but with low accuracy. We also note that the differences in SPT are larger for increasing polymer 
reservoir packing fraction and for size ratio <?=0.6. The SPT performance is worse when the number of polymer in 
the mixture is relatively high compared to the number of colloids. 



In this section we present results on the liquid-gas interfacial tension of AOV colloid-polymer mixtures of size ratio 
5=1, using three independent simulation techniques explained above. First we use the scaling of the interfacial width. 
The Fig. shows typical snapshots of the liquid-gas interface for rf v =0.95, 1.05, 1.4, and 2.0. One observes clearly 
that the difference in densities of the two coexisting phases increases with increasing rf v . Moreover, the liquid-gas 
interface becomes sharper upon increasing jf v . The interfacial order parameter ^f(z) is determined and fitted to an 
error function to extract the interfacial width. In Fig. we plot the square of the interfacial width as a function of the 
logarithm of the lateral dimension L. Employing a linear fit to our data and using Eq. (|44|) . we determine 7i g . In Fig. 
|SJi and Fig. [Hp we plot the liquid-gas interfacial tension as a function of rf c — rf c and rf v respectively. We determine 
the colloid densities rf c and rf c at coexistence from the measured density profile p c (z). 

Than in the second method we determined the probability distributions P{N c )\ Z( , iV t for 77^ =0.9, 1.05, 1.15, 1.5, 
2.0, 3.0, and 4.0 in a cubic simulation box of length L = 7, 8, 9, 10, and 11. We show the probability distribution for 
box length L — 10 for varying rf v in Fig. [5] We observe clearly that the density jump at coexistence and 7i g .L, i.e., 
the difference of the maxima and minimum, increases with rf v . We determined 7i g for an infinite system by plotting 
7i gi i as a function of 1/L 2 in Fig. ITUI and extrapolating, L — + 00, according to Eq. (JSJ). Again, we plot the liquid-gas 
interfacial tension as a function of if c — r\ 9 c and r\ r v in Figs. and Fig. [5p, respectively. We determined if c and rf c at 
coexistence from 



where the factor 2 arises from the normalization l|47|l . Finally, we also employed Young's equation to obtain an 
estimate for 7i g . Dijkstra et al. 01 studied the wetting behavior of AOV model colloid-polymer mixtures using 
computer simulations. From measuring adsorption isotherms they concluded that the colloidal liquid phase wets the 
wall completely upon approaching the gas branch of the gas- liquid bulk binodal for values of rf v < 1.05. In this 
region the contact angle vanishes and 7i g = (7 wg — Twi), where 7 wg and 7 w i are wall-fluid tensions 7 w f computed at 
liquid-gas coexistence. We hence carried out simulations of the wall-fluid tension for statepoints at coexistence below 
the wetting transition point: rf p = 0.935, and 0.977 and also slightly above the wetting transition point: rf v — 1.14 
and 1.25. The coexisting densities were previously determined in Gibbs ensemble simulations (see Fig. |3J). We stress 
that the results obtained on the basis of Young's equation should be taken with great care, as the wall-gas tension in 
our simulations were obtained from a gas phase in contact with a planar hard wall, which is a metastable state with 
respect to the equilibrium state of a macroscopic wetting layer adsorbed at the gas-wall interface. 

However, DFT results show that the contact angle differs only less than few a percent when the metastable gas- wall 
density profile is employed instead of the stable density profile that includes the wetting layer. In Fig. [S^i. and Fig.|Sp 
we plot the difference 7 wg — 7 w i for rf v = 0.935, 0.977, 1.14, and 1.25. Comparing the results obtained for 7i g from the 
three different routes, we find good agreement. The agreement with DFT results is good for high values of rf c — rf c 
and rjp, but deviates close to the critical point, as might be expected. Similar deviations were found for Tig, by Vink 
et al. using the probability distribution method in simulation of the full mixture for size ratio q = 0.8 |22|. 



C. Liquid-gas interfacial tension 




(53) 
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V. CONCLUSIONS 

In conclusion we investigated the wall-fluid tension of the AOV model colloid-polymer mixtures of size ratio q = 
0.6 and q = 1 using Monte Carlo computer simulations. We used a thermodynamic integration method and a 
shoulder potential approximation for the hard-core potentials to determine the free energy of the bulk system and the 
inhomogeneous system. The wall-fluid interfacial tension is the surface excess free energy per unit area, and is in good 
agreement with the DFT results. The SPT wall-fluid interfacial tension is in overall agreement with simulations, but 
the comparison is worse for increasing polymer reservoir packing fraction. We also investigated the colloid and polymer 
adsorption of the colloid-polymer mixture at a planar hard wall and we found good agreement with DFT results. We 
derived a SPT expression for the adsorption of colloid polymer mixtures at a hard wall. The expression reproduce the 
essential features of the adsorption, but with low accuracy. In addition, we studied the liquid-gas interfacial tension 
of the AOV model colloid-polymer mixtures of size ratio q = 1 using i) the dependence of the interfacial width on 
the logarithm of the lateral size of the simulation box as predicted by the capillary wave theory, ii) the probability 
distribution method, and iii) Young's equation in the complete wetting regime. We find remarkably good agreement 
between the different sets of results. Moreover, as we used the effective one-component simulations in the probability 
distribution method, we were able to investigate the interfacial tension for high rf v . As the bulk binodal from DFT, 
equivalent to that of free- volume theory [ll|, agrees well with the simulation results for rf v > 1.5, one might expect a 
similar agreement for the interfacial tension. We found that our simulation results for 7i g approaches the DFT results 
upon increasing 77^ and agree well with the DFT results only for 77^ > 3. Close to the critical point, deviations are 
found between the simulation and DFT results as expected due to the shift in critical point. 

The good agreement of our simulation results for wall-fluid and fluid-fluid interfacial tensions with those of the DFT 
hint at a reliable description of the magnitude of contact angle of the colloidal liquid-gas interface and a hard wall 
[l3|. Nevertheless, carrying out detailed simulations at bulk coexistence in the (numerically very demanding) partial 
wetting regime of large rf p is an interesting issue for future work, in particular in light of the fact that the contact 
angle can be measured in a direct (though not easy) way in experiments 0, ITcj . Furthermore our thermodynamic 
integration technique is well suited to determine the free energy of confined crystals and hence to predict the full 
phase behavior of confined colloid-polymer mixtures; work along these lines is in progress. Finally we mention that 
interfacial properties at curved substrates have attracted recent interest 0, |E(J ; colloid-polymer mixtures are well- 
suited to investigate such situations [To| . 
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FIG. 1: The average number of overlaps ( — ) A , between colloidal particles and the walls per unit volume as a function of the 
height of the step potential for the pure colloidal system (r/p — 0) with packing fraction r\ c = 0.4. For f3\ > 10 the number of 
overlaps is zero within statistical fluctuations and the system behaves like the hard-sphere system. The line is a guide to the 
eye. 




•c 

FIG. 2: The reduced wall-fluid interfacial tension P'ya 2 of hard spheres adsorbed at a hard wall as a function of the colloidal 
packing fraction r/ c . We compare our simulation results (open circles) with Monte Carlo simulations |27| (open squares) and 
Molecular Dynamics simulations |2q | (crosses). The dotted line indicates the result from SPT and the solid line denotes the 
DFT result. 
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FIG. 3: Phase diagram of the AOV model for size ratio q = 1 as obtained from simulations, taken from Refs. |l7| (symbols), 
and free volume theory |Tl| (dashed line) as a function of the colloid packing fraction r\ c and the polymer reservoir packing 
fraction r[ v . F and S denote the stable fluid and solid (fee) phase. F+S and F+F denote, respectively, the stable fluid-fluid, 
and fluid-solid coexistence region. 
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FIG. 4: The wall-fluid interfacial tension of the model colloid-polymer mixture adsorbed against a hard wall. The symbols 
denote simulation results, dotted curves denote SPT results 0], the solid curves denote DFT results Q|. a) Size ratio q = 0.6 
and rjp =0, 0.2, and 0.4; b) size ratio q = 0.6 and rjp =0.5, 0.6, and 0.7; c) size ratio q — 1 and rf v =0, 0.2, and 0.4; d) size ratio 
q = 1 and rf v =0.7, 0.9, and 1.0. 



18 




0.05 



-0.05 - 



-0.1 - 




0.4 0.5 




0.4 0.5 



FIG. 5: The adsorption Fa^ of the colloid-polymer mixture at a hard wall as a function of the average colloid packing fraction 
(rjc)- Simulation results for polymer reservoir packing fraction rf v = (open triangles), rf v = 0.2 (open circles), and rf v = 0.4 
(open squares) are compared with results from DFT (solid lines) and SPT (dashed lines). The DFT results for rf p = are 
omitted for clarity, a) Colloid adsorption for size ratio g=0.6; b) Polymer adsorption for q=0.6; c) Colloid adsorption for q=l; 
d) Polymer adsorption for q=l. 
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FIG. 6: Typical snapshots of the colloid coordinates, obtained from Monte-Carlo simulations based on the exact effective 
one-component Hamiltonian for the colloids, of the liquid-gas interface of the model colloid-polymer mixture with size ratio 
q = 1 and a) polymer reservoir packing fraction rf v — 0.95; b) rf v — 1.05; c) rf v = 1.4; d) rf v — 2. 
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FIG. 7: The squared width of the liquid-gas interface of the model colloid-polymer mixture with size ratio q — 1 for polymer 
reservoir packing fractions rf v =0.95, 1.05, 1.4, 2 (from top to bottom), as a function of the logarithm of the lateral dimension 
of the interface. Symbols denote simulation results, lines indicate linear fits to the data. 
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FIG. 8: The reduced liquid-gas interfacial tension f3^<J^ of the model colloid-polymer mixtures with size ratio q = 1. Open 
squares denote simulation results from the capillary wave method. Filled squares indicate simulation results using the probability 
distribution method. Open circles refer to simulation results employing the difference 7 wg — 7 w i. The solid curves indicates the 
DFT results. The star indicate the position of the critical point, rf v = 0.86, obtained from simulations, a) As a function of the 
difference in gas and liquid packing fractions r\ c — jf c at coexistence, b) As a function of the polymer reservoir packing fraction 
r/p, at coexistence. The inset shows the data close to the critical point on an enlarged scale. 




FIG. 9: Logarithm of the probability P{N C ) (not normalized, as P(0) is taken to be 1) of finding N c colloids with diameter 
<7 C for a colloid-polymer mixture with size ratio q = 1 in a cubic box with volume V = lOOOo-j? at varying polymer reservoir 
packing fraction rf v =0.9, 1.05, 1.15, 1.5, and 2 (from top to bottom), as a function of r) c . All statepoints are at coexistence. 
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FIG. 10: The system size dependence of the (scaled) liquid-gas interfacial tension /?7i g , l<7c of a colloid-polymer mixture with size 
ratio q = a p /a c — 1 in a cubic box with volume V = L 3 al at varying polymer reservoir packing fraction rf v = 0.9, 1.05, 1.15, 1.5, 
and 2.0 (from bottom to top). The dashed lines are least-squares linear fits. 



